Dynamical localization, measurements and quantum computing 
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We study numerically the effects of measurements on dynamical localization in the kicked rotator 
model simulated on a quantum computer. Contrary to the previous studies, which showed that 
measurements induce a diffusive probability spreading, our results demonstrate that localization 
can be preserved for repeated single-qubit measurements. We detect a transition from a localized to 
a delocalized phase, depending on the system parameters and on the choice of the measured qubit. 
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In 1979 the dynamical localization of quantum chaos 
was discovered in numerical simulations of the kicked ro- 
tator model . It was found that the unbounded classi- 
cal diffusion typical of chaotic dynamics is suppressed by 
quantum interference effects This interesting phe- 

nomenon found its explanation on the basis of an analogy 
with the Anderson localization in disordered lattices |3( 
(see also |4|). Manifestations of dynamical localization 
appear in various physical systems. Its first experimen- 
tal observation was obtained with hydrogen and Ryd- 
berg atoms in a microwave field 0. Recently, a signifi- 
cant technological progress in manipulating cold atoms 
by laser fields allowed to build up experimentally the 
kicked rotator model and to study dynamical localiza- 
tion in real systems in great detail [E 0> HI • 

Since localization appears due to quantum interference 
it is natural to expect that it is rather fragile and sensitive 
to noise and interactions with the environment. Indeed, 
in theoretical and experimental studies it was shown that 
even a small amount of noise destroys coherence and lo- 
calization 0, H, • Measurements represent another type 
of coupling to the environment 0, and it is of funda- 
mental importance to understand their effects on dynam- 
ical localization. Theoretical and numerical studies show 
that measurements destroy localization and induce a dif- 
fusive energy gro wth like in the case of a noisy environ- 
ment [ill 1 1 2L I X 3l - For weak continuous measurements, 
discussed in the rate of this growth can be much 
smaller than the diffusion rate induced by classical chaos. 
However in the limit of strong coupling to the measure- 
ment device the quantum diffusion rate becomes close to 
its classical value. A similar situation takes place in the 
case of projective measurements, considered in HI EH. 

The interest in measurement procedures enormously 
grew up in the last years due to progress in quantum in- 
formation processing Indeed, the extraction of infor- 
mation from a quantum computation is always reduced 
to a final measurement of the quantum register. Vari- 
ous experimental implementations were discussed for the 
realization of the readout procedure in quantum optics 
systems nam and solid state devices 1171 
Moreover it has been shown that a quantum computation 



can be performed completely by a sequence of measure- 
ments applied to an initially entangled state [ilj . At the 
same time, measurements represent an important part of 
various quantum algorithms, including the famous Shor 
algorithm for the factorization of integers • Therefore 
it is important to investigate the effects of measurements 
on quantum computers operating nontrivial algorithms. 

An interesting example is the quantum algorithm pro- 
posed in which allows to simulate the evolution of the 
kicked rotator on a quantum computer. This algorithm 
essentially uses the Quantum Fourier Transform (QFT) 
and controlled phase gates. It realizes one map itera- 
tion in a polynomial number of quantum gates (0(n^)) 
for a wave vector of size N — 2 n " . Here n q is the num- 
ber of qubits (two-level quantum systems) onto which a 
kicked rotator wave function is encoded. For moderate 
kick amplitudes, this algorithm can be replaced by an 
approximate one which uses all the qubits in an optimal 
way and performs one map iteration in 0(n^) elemen- 
tary gates 24] . In this form the algorithm can simulate 
complex dynamics, e.g. the Anderson transition, with 
only a few (~ 7) qubits. This makes it accessible for 
possible future realization on NMR based quantum com- 
puters. Indeed, all the elements of the algorithm have 
already been implemented on NMR quantum computers 
[25I l26| . Therefore it represents an interesting testing 
ground for the investigation of the measurement effects 
on dynamical localization in a quantum computation. 

The quantum evolution of the kicked rotator is de- 
scribed by the unitary operator U acting on the wave 
function ip |4|: 



(1) 



Here yj is the wave function after one map iteration, Uk 
represents the effects of the kick in the phase represen- 
tation and Ut describes the free rotation in the momen- 
tum basis n with h = —id/dO (we use units with h = 1). 
The dimensionless parameters fc, T determine the kick 
strength and the rotation phases, so that the classical 
limit corresponds to k — » 00, T — > with the chaos pa- 
rameter K = kT constant. Here we study the regime of 
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dynamical localization corresponding to I <C N, where 
/ « k /2 is the localization length |^|. 

The quantum algorithm simulating the evolution ifTJl 
operates as described in [2^, 0|- The wave function 

in the momentum representation with N — 2™« lev- 
els is encoded on a quantum computer with n q qubits. 
In this way n = —N/2 + j where the index j — 
0, . . . , N — 1 is written in the binary representation as 
j = (ai,a 2 ,...,a m ,...,a nq ), with a m = or 1. As 
the initial state we choose the momentum cigenstate at 
n = uq — 0, which can be efficiently prepared from 
the ground state. Then, as described in |23j the rota- 
tion Ut is performed in 0(n 2 q ) controlled-phase gates. 
After that the QFT transforms the wave function to 8- 
representation in 0(n 2 ) quantum elementary gates (see 

(l4|). The kick operator Uf. is realized in 0(n q ) gates 
with the help of an additional register j2^| , or, for mod- 
erate k values, it can be approximately implemented in 
0(n q ) gates without any ancilla, following |2J|. Finally, 
ip is transformed back to the momentum basis by the 
inverse QFT. Here we assume that the gates are imple- 
mented without errors, keeping the analysis of imperfec- 
tion effects for further studies. 

To study the effects of measurements on the dynam- 
ics given by the above algorithm we assume that after 
each map iteration Q a projective measurement of a 
chosen qubit m is performed. The measurement can 
be represented as the action of two projection opera- 
tors Po(m) and Pi (to) giving for a m an outcome or 

1 with the probability ||Po(m)-0|| 2 or ||Pi(to)V>|| 2 . re- 
spectively. The measurement projects the wave func- 
tion onto one of two subspaces of the total Hilbert space, 
corresponding to momentum states labeled by the in- 
dexes j = (ai, 02, ■ ■ • , a m , ■ ■ ■ , a n ) with fixed a m = or 
a m = 1. Each subspace is composed of N/2 states, given 
by the direct sum of 2 m_1 cells of L = 2 nq ~ m consec- 
utive momentum states. For example, for m = 1 the 
most significant qubit is measured and tp is projected 
onto momentum states with —N/2 < n < (a% = 0) or 
< n < N/2 — 1 (ai = 1); for m — n q the least signifi- 
cant qubit is measured and ip is projected onto even and 
odd momentum states. Such a measurement is the most 
natural one for the quantum computation process. 

Thus, the evolution with measurements is given by the 
following equation for the density matrix p 

fj = P (m) -U k -U T -p-U T - Ut ■ Pa(m) + 

Pi (to) • U k ■ U T ■ p ■ U T ■ Ut • Pi(m) (2) 

Here, p is the density matrix after one map iteration with 
measurement. The direct simulation of this equation is 
quite costly, since N 2 components should be iterated. To 
avoid this difficulty we used the method of quantum tra- 
jectories [i^ • In this method for one quantum trajectory 
the wave function ip evolves according to (JIJ; after each 
map iteration ip in the momentum representation is pro- 
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FIG. 1: (Color on line) Dependence of the second moment 
(n 2 ) on time t. Here T — 2, k — 2 and color marks the 
value of n q . The upper group of curves corresponds to the 
measurements of the least significant qubit m — n q for n q = 
12 (blue), 11 (violet), 10 (red) and 9 (black), from top to 
bottom. In the lower group one of most significant qubits 
is measured with m = n q — 8 for n q = 12 (black) and 9 
(yellow) (data for n q — 10, 11 give same superimposed curves 
and are not shown). The lowest green fluctuating curve is the 
evolution without measurements. The dashed line shows the 
diffusive growth (n 2 ) ~ t. The inset shows the dependence of 
the IPR (oat (colors are as in the main plot, the dashed line 
shows the diffusive growth £ ~ Vi). 



jected on the subspaces with a m = or 1 according to 
the probability | \Po(m)ip\ | 2 or |Pi(m)-0|| 2 , respectively. 
After the renormalization, this gives the wave function 
ip n in the momentum basis. The density matrix and the 
expectation values of observables are then obtained by 
averaging over M quantum trajectories. 

To characterize the quantum evolution with measure- 
ments we compute the following quantities: the prob- 
ability distribution p nn w (|?/Vt| 2 ), obtained by averag- 
ing \ip n \ 2 over M quantum trajectories; the second mo- 
ment of the probability distribution, given by (h 2 ) — 
Tr{h 2 p) ^ n n 2 (\ip n \ 2 ) ; the Inverse Participation Ra- 
tio (IPR) £ = 1/ p 2 nn « 1/ £„ |(|V„| 2 )| 2 which deter- 
mines the number of states on which the average proba- 
bility is distributed. Within statistical fluctuations these 
quantities remain unchanged for a variation of M from 
20 to 500 and we represent them for M = 50. 

The dependence of (h 2 ) and £ on the number of map 
iterations t is displayed in Fig^for different n q and to. 
The probability distribution (|"0n| 2 ) for n q — 10 is shown 
in Fig[21 These data clearly show that the measurement 
of the least significant qubit completely destroys localiza- 
tion generating a diffusive behavior. Indeed, the second 
moment (n 2 ) and the IPR £ grow diffusively (see FigQ 
up to spread of the probability over the whole computa- 
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FIG. 2: (Color on line) Probability distribution for k = 2, T = 
2,n q = 10 at t — 5 x 10 J . Measurements are done for 
m — n q — 8 that preserves localization (black curve) and for 
m = n q that leads to extended distribution (red curve). The 
distribution for evolution without measurements is shown by 
green/gray curve. Data are averaged over 50 quantum trajec- 
tories. The inset shows \^ n \ 2 for a single quantum trajectory 
(same colors). 
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FIG. 3: Dependence of the averaged IPR (£) on k for mea- 
surements of one of the most significant qubits m = n q — 8, 
for n q — 9 (squares); 10 (diamonds); 11 (triangles); 12 (stars). 
IPR values are averaged over 1000 kicks around t = 5 x 10 5 ; 
T = 2. The inset shows the same dependence in the absence 
of measurements. 



tional basis, as shown in Fig|21 The extended distribu- 
tion (l^nl 2 ) is formed by a superposition of probabilities 
\ip n \ 2 generated by single quantum trajectories (see Fig. 
13 inset, which shows that each \ip n \ 2 is relatively nar- 
row). On the contrary, the measurement of one of the 
most significant qubits does not destroy localization, as 
clearly illustrated in Figs 121 This stri king result is very 
different from the previous studies ll'2t Il3| where lo- 
calization was always destroyed by measurements. 

To understand the origin of this behavior we investi- 
gate the dependence of the averaged IPR (£) on the kick 
amplitude k for different number of qubits n q (see Figl^Jl. 
For k < k c ss 6 the IPR is independent of n q , correspond- 
ing to a localized regime. On the contrary, for k > k c the 
IPR starts to grow with the system size N = 2 7lq , indi- 
cating a transition to a delocalized phase. We explain the 
appearance of this transition in the following way. The 
measurement process determines the cell size L = 2 nq ~ m 
inside which the coherence of quantum dynamics is pre- 
served. If the unperturbed localization length I is much 
smaller than the cell size, then measurements do not de- 
stroy dynamical localization. While, if I ^> L, the wave 
function propagates over different cells, measurements 
destroys quantum coherence between nearby cells and 
this leads to a diffusive propagation over the computa- 
tional basis. 

According to our data the delocalization transition 
takes place when 

£o«2Z«fc 2 «£/5 (3) 

where £o is the IPR for the dynamics without mea- 
surements (see the inset in FigEJ). This relation shows 
that the transition can be obtained by tuning k at fixed 
n q — m or by an appropriate variation of m at fixed k. 



Our numerical data confirm this estimate. Indeed, for 
m = Uq — 9 we obtain that k = 10 is localized, while at 
k = 12 delocalization takes place (data not shown). It 
is interesting to note that the oscillations of (£) in Fig|3| 
are correlated with the oscillations of £o> thus confirm- 
ing that the delocalization border is determined by the 
unperturbed localization length I (these oscillations are 
produced by dynamical correlations which affect the clas- 
sical/quantum diffusion rate related to the localization 
length I as discussed in |5g]). 

To study the quantum dynamics at larger time scales 
we use the random quantum phase method proposed in 
p^ . It is based on the fact that after a projection on 
a given quantum state induced by a measurement the 
quantum phase is not defined. Therefore one can as- 
sume that states associated to different outcomes of the 
measurement procedure have a random relative quantum 
phase. Thus, after a measurement of the m-th qubit 
the state \<j>) is replaced by e ifi ° Pa{m)\4>) + e if3l Pi(m)|(/>), 
where the phases /?o,i are random. This approach allows 
to reduce significantly the computational cost of the sim- 
ulation, since it effectively integrates the dynamics over 
many quantum trajectories. 

The comparison of the two computational methods is 
presented in Fig0] for diffusive, localized and critical 
regimes. Both methods give consistent results for (n 2 ) 
(Fig0J and the IPR (data not shown). With the ran- 
dom quantum phase method we can follow the evolution 
for very large times (up to t = 10 7 ) at which localiza- 
tion is still preserved (see FigQ]b). This computational 
method allows also to understand in a better way why 
localization is not destroyed by measurements. Indeed, 
the effects of random phase fluctuations /?o,i appear only 
at the cell boundaries. Hence, for L 3> I they do not 
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FIG. 4: (Color on line) Time dependence of the second mo- 
ment (n 2 ) of the quantum distribution obtained by the com- 
putation with quantum trajectories (blue and black curves) 
and with the random quantum phase method (yellow/gray 
and green/gray curves); T — 2, n q = 10. Panel a) shows 
diffusive regime for k — 6 (upper blue and yellow curves) and 
k = 2 (lower black and green curves), for m — n q ; the straight 
line gives the diffusive law {n 2 } ~ t. Panel b) shows a local- 
ized regime for k = 2 (lower curves) and a near critical case 
for k — 6 (upper curves), for m = 2; the straight line shows 
anomalous diffusion (n 2 ) ~ t ' 2 , colors are as in panel a). 



affect the momentum states located on a distance larger 
than I from edges and localization is preserved j2t|. We 
think that the same mechanism qualitatively explains the 
results obtained in [^Ol where it was found that measure- 
ments of a 1 /2-spin detector coupled to the kicked rotator 
do not destroy localization. In this case the effective cell 
size L is the total number of rotator momentum states 
and thus localization is preserved since / -C L. 

In conclusion, we studied the effects of measurements 
on dynamical localization in a quantum algorithm sim- 
ulating the kicked rotator. Contrary to the common 
lore the localization is not always destroyed by measure- 
ments, and a transition from localized to diffusive dy- 
namics takes place when system parameters are varied. 
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